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Abstract 

We study the asymptotic shape of self-avoiding random walks (SAW) on the 
backbone of the incipient percolation cluster in d-dimensional lattices analyt- 
ically. It is generally accepted that the configurational averaged probability 
distribution function {P&{r,N)) for the end-to-end distance r of an N step 
SAW behaves as a power law for r — > 0. In this work, we determine the cor- 
responding exponent using scaling arguments, and show that our suggested 
'generalized des Cloizeaux' expression for the exponent is in excellent agree- 
ment with exact enumeration results in two and three dimensions. 



PACS: 5.40.-a, 61.41. +e, 61.43.-j 
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Linear polymers have a broad spectrum of applications in such different fields as biology, 
chemistry, and physics, see e.g. and theoretical work has been extensively documented 

in the literature More recently, much interest has been drawn to the effects that 

environmental disorder may have on the configurational properties of such linear structures 
|7]-|9|]. As a paradigmatic model of linear polymers, self-avoiding random walks (SAW) have 
captured the attention of researchers for decades and are currently extensively studied J| ||. 
Structural disorder is generally modelled by percolation [TO,0], and interesting effects are 
expected at the percolation threshold || . Since asymptotically long SAW can only exist on 
the backbone of the cluster, in which all dangling ends have been eliminated, one may study 
SAW directly on the backbone. For relatively short chains, one may still expect that the 
statistical properties of SAW do not depend on whether they are studied on the backbone 
alone or on the corresponding percolation cluster. In what follows, we adopt this point of 
view and study SAW directly on the backbone of the cluster. 

On percolation clusters, one can distinguish between two different metrics, the usual 
Pythagorian one, or r-space metric, and the topological one, or £-space metric, which have 
been extensively discussed in the literature p}fL3|j. The corresponding length scales are 
related by i ~ r dmin , where d m i n is the fractal dimension of the shortest path [11]. 



The scaling behavior of SAW is characterized by the mean end-to-end distance after 
N steps. On structurally disordered systems such as percolation clusters, two types of 
averages need to be performed. One first averages over all SAW configurations of N steps 
on a single backbone starting from the same origin, as obtained e.g. by exact enumeration 
techniques fl4 |, followed by a configurational average over many backbones. In r-space, the 



former average is denoted as r 2 (N), and the second one by (r 2 (N)) ~ jV 2 ^, which defines 
the exponent v r . Correspondingly, in £-space one has, (i(N)) ~ N Vt , where v t = v r G? m i n . 

More generally, one is interested in the configurational averaged probability distribution 
function (PDF) for the end-to-end distance on the backbone. In r-space, it is denoted as 
(Pb(^, N)), and is expected to obey the scaling form 

(PB(r,N)) = -L w f r (r/N^) for N » 1 , (1) 

which is normalized according to dr r d ? _1 (-Pb^ N)) = 1, where df is the fractal dimen- 
sion of the backbone. The scaling function f r (x) behaves as f r {x) ~ x 9 ^, for i<1, and in 
this limit Eq. (|1]) becomes 

(ft(^))~ifi) 91 ^ r«i\^ , (2) 



where g[ is the scaling exponent to which we draw our attention here. In £- 
space, the PDF obeys (Pb(£,N)) = N~ Ued ? fi(£/N Ut ) and is normalized according to 
/ °° d£ {Pb(£, N)) = 1, with df being the fractal dimension of the backbone in £-space, 
df = df/d m i n . The corresponding exponent in the case £ <C N Vt is denoted by g[ and is 
related to g[ by g{ = g[/d m i n (see e.g. [ jl5l , fjj| ) ■ The interest in studying the PDF (Pb(£, N)) 



is because fluctuations are definitely smaller in £-space than in r-space, and more accurate 
results for g\ can be derived by determining g[ and afterwards using the relation g\ = g\d m i n . 
In the case of regular lattices (g[ = g\ = gi), des Cloizeaux [|I7[] showed that g\ = 



(7 — l)/z/p, where z/p = 3/(d + 2) is the Flory exponent |18| and 7 describes the total number 
Cn of SAW configurations of steps, 
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C N = fi N N^ 1 for N > 1 , (3) 

where \x is the effective coordination number of the lattice (see e.g. ||). Numerical results 
for SAW on the backbone of percolation clusters at criticality in two dimensions seem to 
be roughly consistent with the form g\ = (7! — V)jv r ||15|| , where 71 is the corresponding 
enhancement exponent on the backbone. More recent results in three dimensions clearly 
disprove this relation, so the question remains whether a generalization of the des Cloizeaux 
expression can still provide an accurate framework for estimating the exponent g\ analyti- 
cally. The goal of the present paper is to provide such a generalization of the des Cloizeaux 
relation valid in all dimensions. 

To this end, we follow de Gennes [|J and study the probability (Pb(o, N)) that the iV-th 
step of the SAW is on a backbone site located at a lattice distance a = 1 from its starting 
point at r = 0. This probability can be written as 

W«=1,^)~ <C "* ( ''° 1)> (4) 

where (Cjv,b(o = 1)) represents the (configurational) average number of SAW that, after N 
steps, arrive at a distance a = 1 from the origin, and (Cn,b) denotes the total number of 
SAW of iV steps. The latter is the configurational average analog of Eq. ([3|), with /x and 7 
replaced by the exponents fi\ and 71. The index 1 reminds us of the underlying multifractal 



behaviour of these two exponents for percolation systems [jlq] (see also [19]). We assume 
the former to be 



(C Nfi (a=l)) = tf [^) ' F B (5) 

where R = (r 2 (A^)) 1 / 2 ~ N Ur , and the factor F-q ~ R~ K < 1 represents the additional 
difficulty of the N step SAW to return close to its starting point due to topological constraints 
caused by the statistical nature of the embedding structure. For regular lattices, this factor 
reduces to unity since in those cases there is no distinction between cluster and lattice 
sites, and the exponent df is replaced by the spatial dimensionality d of the lattice ||. For 
deterministic fractals such as Sierpinski gaskets, the factor Fb is also expected to be unity, 
since in those cases the fractal structure embedding the SAW is also unique j2U| . 

The exponent k can be determined by considering first the probability that an arbitrary 
site belongs to the infinite cluster, P^. Near the percolation threshold p c , P^ behaves as a 
function of the concentration p (> p c ) of occupied sites as [[njjnj 

Poo~(p-Pe) P ~C t,/ " (6) 

where the correlation length £, diverging as £ ~ \p — p c \~ v near p c , is a measure of the linear 
size of the finite clusters in the system (or similarly, the mean distance between two sites 
on the same finite cluster), and v is the correlation length exponent. Eq. (|6|) suggests that 
the probability P 00 (L), to find a site belonging to the incipient cluster within a distance L 
from a given cluster site (and consequently at least one path connecting these two sites), 
scales as P^L) ~ L~PI V for L < £ |11| . Identifying the length scale L with the mean size R 
of the SAW, L ~ R, we obtain Fb ~ R~^^ h ', i.e. k = P/u. Accepted values for the critical 
exponents reported so far ]p!0| - [T3|j2l| -p5| are summarized in Table |. 
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Using Eqs. (||), (||) and the above result for Fb, we find 



^a = ^N))~ - f ^ ri ~^{ Wr ) , (7) 



and comparison with Eq. (|2|) yields 



9 ; = ^ + ^ (8) 



which is denoted as the generalized des Cloizeaux relation. 

Numerical calculations based on exact enumeration techniques in two and three dimen- 
sions (taken from Ref. |16[]) are reported in Fig. [I] for (P^(£,N)) |[26|| . The corresponding 
numerical values for g\ and g\ are summarized in Table f|. The suggested relation Eq. (§) 
is in excellent agreement with the numerical data. Additionally, it should be noted that 
Eq. (H) yields the correct value expected when d > 6, i.e. g[ = 2. Indeed, when d > 6 the 
backbone of the incipient percolation cluster becomes topologically one dimensional, since 
loops are irrelevant on large length scales, and one simply has 71 = 1. Thus, g[ is solely 
determined by the second term in Eq. (p. The value g\ = 2 for d > 6 can be understood 
by noting that in this case Fb ~ (a/R) df with df = 2 for d > 6, since SAW are equivalent 
to the embedding backbone and the probability to return close to the starting point just 
decreases as R~ 2 , i.e. g\ = 2. The good agreement between Eq. (§) and the numerical 
results supports the ansatz made in Eq. (^). 

Financial support from the German-Israeli Foundation (GIF), the Minerva Center for 
Mesoscopics, Fractals, and Neural Networks and the Deutsche Forschungsgemeinschaft is 
gratefully acknowledged. 
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FIGURES 



FIG. 1. Scaling plots of the PDF in ^-space on the backbone of the incipient percolation cluster 
in (a),(b) d = 2 and (c),(d) d = 3, from the exact enumeration results of Ref. fl6[j . In (a) and 
(c), the whole PDF is shown, plotted as £ d ® (P B {£,N)) versus t/N" 1 , for N = 39 (full circles) 
and 40 (full squares), averaged over 5 ■ 10 3 backbone configurations. In (b) and (d), the region 
with £ < N"* is shown in more detail, plotted as £ d e (P B (£,N)) /{£/N Ui ) d * = (N^ f* {P B (£,N)) 
versus £/N Ue 26]. Case d = 2: In (a) the continuous line corresponds to the case £ 3> N" e and is 
not further discussed here (see [|l6| for details). For £ <C N Ul the dashed line in (b) has the slope 



0.49±0.05 and is a fit with the ansatz (A^) d F (p B (£,N)) ~ {£jN v ^ for £ < N VI (cf. Eq. (|) and 
text), in very good agreement with the prediction g[ = (71 — l)/ve + 0/(^d m i n ) = 0.474 (derived 
from Eq. (||) using the relation g\ = g\/d m i n ). The dotted line is for illustration only and has the 
slope 0.382, resulting from the relation (71 — l)/vi — 0.382, in disagreement with the numerical 
results. Case d = 3: Similarly to (b), in (d) the dashed line has the slope 0.67 ±0.05, in excellent 
agreement with the prediction g[ = (71 — + /?/(fdmin) — 0.666. The dotted line has the slope 
0.32, resulting from (71 — l)/z^ — 0.32, clearly revealing the inadequacy of this relation to fit the 
numerical results. 
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TABLES 





d = 2 


d = 3 


d > 6 


p 


5/36 a 


0.417 ±0.003 b 


l c 


V 


4/3 a 


0.875 ± 0.008 b 


1/2 C 




91/48 a 


2.524 ± 0.008 b 


4 d 


df 


1.6432 ± 0.0008 c 


1.87 ±0.03' 


2 d 


dmin 


1.1306 ±0.0003 g 


1.374 ± 0.004 h 


2 d 


df 


1.446 ± 0.001 1 


1.36 ±0.02' 


l d 



TABLE I. Critical exponents and fractal dimensions for the incipient percolation cluster, for 
spatial dimensions d = 2, 3, and d > 6, from a exact results Ref. [^IJ, b Monte Carlo simulations 
Ref. ||, c exact results Ref. ||, d Refs. 00 > CMonte Carl ° simulations Ref. ||, f Monte Carlo 
simulations Ref. p5| , g Monte Carlo simulations Ref. [13], 
'the relation 



g, 1 
df/d n 



1 Monte Carlo simulations Ref. [12], and 



Here, df = d — (3/v is the fractal dimension of percolation clusters. 





d = 2 


d = 3 


d > 6 


v r 


0.787 ±0.010 


0.662 ± 0.006 


1/2 


7i 


1.34 ± 0.05 


1.29 ±0.05 


1 


9l 


0.56 ±0.10 


0.90 ±0.10 


2 


A 


0.49 ± 0.05 


0.67 ±0.05 


1 


9l = 9l d m in 


0.55 ± 0.06 


0.92 ± 0.08 


2 


g[ (present conjecture) 


0.54 ± 0.07 


0.916 ±0.080 


2 



TABLE II. Results for SAW on the backbone of critical percolation clusters (from Ref. fllGf). 
Note that the value of g[ slightly differs from the value given in Ref. [p3[| , as we determined it 
more accurately here by plotting (Pq(£,N)) as shown in Fig. p]. The last line displays the present 
theoretical results from the generalized des Cloizeaux relation g\ = (71 — l)/u r + (3/v, Eq. (^). 
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